Suppression of current in transport through parallel double quantum dots 



Tae-Suk Kim a ' b and S. Hershfield 6 
a Institute of Physics and Applied Physics, Yonsei University, Seoul 120-749, Korea 
b Department of Physics, University of Florida, Gainesville FL 32611-8440 
^ [ (February 1, 2008) 

We report our study of the I-V curves in the transport through the quantum dot when an 
additional quantum dot lying in the Kondo regime is side-connected to it. Due to the Kondo 
scattering off the effective spin on a side-connected quantum dot the conductance is suppressed at 
low temperatures and at low source-drain bias voltages. This zero-bias anomaly is understood as 
, , ■ enhanced Kondo scattering with decreasing temperature. 
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1 !• INTRODUCTION 

Advances in nanotcchnology have made it possible to fabricate quantum dots or artificial atoms and to study their 
transport properties. Due to the confinement and the Coulomb interaction in the quantum dot, charge and spin can 
be quantized. Quantization of electron charge leads to the almost periodic variation of conductance as a function of 
the number of electrons N or the gate voltage. A conductance peak is observed whenever the charge fluctuation is 
allowed or the N- and (N + l)-states become energetically degenerate. Otherwise electrons in the electrodes cannot 
hop into a quantum dot due to the strong Coulomb repulsion. This Coulomb blockade leads to vanishing conductance 
between two neighboring conductance peaks. Novel Kondo effects can enhance the conductance in the Coulomb 
blockade region when the number of electrons at the quantum dot is odd and at least one electron spin is unpaired. 
The unpaired 5=1/2 in the quantum dot caa be screened by the electrons in the external electrodes connected to 
the quantum dot. This theoretical prgdictionEl of the novel nonequilibrium Kondo effects in quantum dot systems 
was confirmed experimentally recentlyatl In the conventional bulk Kondo systems where magnetic ions are doped 
into normal metals the study of the Kondo effects was confined to the equilibrium state. In quantum dot systems the 
nonequilibrium situation is easily controlled and model parameters can be varied by adjusting the gate voltages. 
. The study of the Kondo effect in the quantum dot systems has been diversified to see the different aspects. The 
' quantum dots are inserted in the Aharonov-Bohm ring to study the effect of many-body interaction on the persistent 
currentO. When electrons are scattered off the magnetic impurity, they experience a phase shift tt/2 at the Fermi 
energy due to the Kondo effect. The direct measurement of this phase shift is not possible in the bulk systems. Using 
a quantum dot inserted in the .Aharonov-Bohm ring the phase shift due to the Kondo effect was studied theoretically!! 
and measured experimentally!! Serial double quantum dots were also studied in the Kondo regimeti Enhanced 
conductance was observeda in an otherwise Coulomb blockade region at the spin singlet-to-triplet transition point 
when the number of electrons on a quantum dot is even. The ground state can be either spin-singlet or spin-triplet 
depending on the geometry of the confinement potential. Perpendicular magnetic fields can induce a transition of the 
ground state between spin-singlet and spin-triplet. This-novcl Kondo effects have been interpreted as coming from 
'"^J ■ the enhanced Kondo temperature at the transition pointed. 

In this paper we study theoretically the transport properties of parallel double quantum dots as depicted in Fig^j. 
One active quantun dot is connected to the source and drain electrodes and the other quantum dot is side- connected 
to the active quantum dot. Since the energy levels in two quantum dots can be controlled separately using the gate 
voltage, different transport regimes can be probed. In this work we present our study of the I-V characteristics of 
this system when the active quantum dot lies in the conductance peak region while the side-connected quantum dot 
in the Coulomb blockade region. In this case charge fluctuations are suppressed and an effective spin 5=1/2 arises 
in a side-connected dot. Electrons passing through the active dot experience the Kondo scattering off an effective 
spin in a side-connected dot. We find that the linear response conductance is suppressed at low temperature due to 
the enhanced Kondo scattering off spin S = 1/2. The spectral function at a side-connected dot develops a Kondo 
resonance peak as the temperature is lowered. On the other hand the spectral function at the active dot becomes 
depleted near lu = with decreasing temperature. We obtained these results using the nonequilibrium Green's 
function method combining with the non-crossing approximation (NCA). 

This paper is organized as follows. In section II we introduce our model Hamiltonian and present the formulation of 
current passing through an active dot. The relevant Green's functions in nonequilibrium are formulated in section III 
using the non-crossing approximation (NCA) which is modified for our system. Numerical results, solving the NCA 
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equations self-consistently, are presented in section IV and conclusion is included in section V. 



II. PARALLEL DOUBLE QUANTUM DOTS 
A. Model Hamiltonian 

Consider the parallel double quantum dots shown in Fig. [I]. In this experimental geometry, the transfer of electrons 
from one electrode to the other is realized by hopping on and off an active quantum dot (the site A). The other dot 
(the site S) is side- connected to the active dot by hopping. Assuming that the energy level spacing is large enough 
compared to the level-broadening due to the tunneling into two electrodes, the model Hamiltonian can be written as 



H — H c i + -ffqdot + Hi , (2.1a 
H e \ 



V Ve pcV c r , (2.1b) 

p=L,R £ Q 



ffqdot= Ys Ei 4A« + E U i n H n U> ( 2 - lc ) 

i=A,S a i=A,S 

H i = E E [ v p$) c UJa* + v;(k) d AaCp , a ] + J2 [wd\ a d Sa + w% a d Aa ] . (2.id) 

p=L,R £ Q a 

The electron creation operator of spin a in the electrode p = L, R is c^- . The index p — L(R) denotes the left (right) 

external electrode. Both electrodes are assumed to be described by the Lorentzian density of states (DOS) with the 
energy dispersion e p £ = eg + pl.r, = ±^eV^ is the chemical potential shift due to the bias voltage applied to the 

electrodes. d\ a is the electron creation operator in a quantum dot i = A, S with A(S) labeling a active (side-connected) 
dot, respectively. Ui is the on-site electron-electron Coulomb interaction in a dot i. 

When U a — Us — 0, two dots become resonant levels and the solution can be found analytically (see Appendix 
A). The model becomes nontrivial when one of on-site Coulomb interactions or both are nonzero. In this paper we 
are going to study the case of Ua — and Us ^ in detail. When Us is the largest of all the model parameters 
and — Es S> \W\, charge fluctuations in a side-connected dot can be neglected and the side-connected dot can be 
treated as an effective spin S ^1/2. In this case removing the charge degrees of freedom in a side-connected dot 
(Schrieffcr-Wolf transformationliil) the interaction between two dots can be derived. 

— - l 2|14^| 2 2|W^| 2 

Has = JasSs ■^2 l d\ a -a a pd AI3 , Jas = „ £ ; g + u s + E s ' *~ 2 ' 2 ^ 

Here Ss = ^2 a p \S;a > ^cf a p < S;f3\ represents the spin degrees of freedom in a side-connected dot in the absence 
of charge fluctuations. Electrons, flowing from the left electrode to the right, at the active dot will experience the 
Kondo scattering off an effective spin of a side-connected dot. Appendix B contains the details of the Ua ^ and 
Us ^ case. 



B. Formulation of Current 



The current operator can be defined as a variation of the number operators of electrons per unit time leaving one 
electrode (p — L,R). 



-e(-N p ) = ±[N P ,B\ = Jr£[^(fc) C ^ AQ -ff, 



(2.3) 



The measured current is given by the thermaL-average of the above current operator and can be written in terms of 
the lesser Green's functions out of equilibriumt 2 ! 



-IE 



E [ V p® G Ac P (k;t,t) - V p *(k) G< A (k;t,t) 



(2.4) 



hex 
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The current can be expressed in terms of Green's functions for the conduction electrons and quantum dots using the 
Dyson equations for the mixed Green's functions 



( TC pk a ^ d Aa(t')) 

V p {k) [ dti G cp {k;t,tx) D A {t x J), 
t 



c 



1 



G Acp (k;t,t') = - (Td Aa {t)c^ ka {t') 



V;{k) / dtt D A (t,h) G cp (k;ti,t'), 



ihD t {t,t l ) = <Td la {t)d\ a {t')>, 1 = Ab- 
using analytic continuation to the real-time axisEl, we can find the lesser part of G Acp (k; t, t') 

Gi cp (k;t,t') = V p *(k) J dh [Dfah) G%{k;ti,t') + D r A (t,h) G<(k;ti,t') 

The current can be expressed as 

I p = ~ ie ^ y E\ v p&)\ 3 \% {^HK p H-s;h] + e<hk(,)-dih]} 

ka 

= - je E / % I^MlM^HK p H-s^)] + s<HKH-s-H]}. 

a 

Here E^ is the self-energy for the active dot A and the following equations are used 



^cp(^j t ) 


= \v p (k) 2 G cp (k;t,t'), 




k 

fde T p (e) 


J it hui — e ± iS ' 




= 2T p {u)f p (±io), 


T p (u) 


= JE \ V ^ 2 



(2.5a) 

(2.5b) 
(2.5c) 

(2.6) 



(2.7) 



(2.8a) 

(2.8b) 
(2.8c) 
(2.8d) 



The above equation for the current is quite general and can be used as a starting point in reducing the expression of 
current in specific cases. 

The average number of electrons at a side-connected dot S should remain the same in a steady state. This condition 
is essential to get the right expression for the current flowing from the left to the right electrodes. That is, the net 
current flowing out of the dot S should vanish for 



^ = ^[ r < d Ld Aa >-W <d Aa d Sa >] = -ieJ2[W*Gi s (t,t)-WG< A (t,t)] 

a a 

The current 1$ was expressed in terms of the mixed Green's functions 

iKG AS (t,t') = <Td Aa (t)dl a (t')>, 
ihG SA (t,t') = <Td Sa (t)d Aa (t') > . 

These mixed Green's functions be written as 

G AS (t,t') = W [ dh D A {t,ti)D s {ti,t') = W f dh D A {t,ti)D s {ti,t'), 
G SA (t,t') = W* [ dh D s {t,ti)D A {ti,t') = W* [ dh D s {t,ti)D A (ti,t'). 



(2.9) 



(2.10a) 
(2.10b) 



(2.11a) 
(2.11b) 
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Here Da and D$ are the fully dressed Green's functions of the dots A and S, respectively. 

ihD A (t,t') = <Td Aa (t)d Aa (t')>, 
ihD 3 (t,t) = <Td Sa (t)4 a (t) > . 



(2.12a) 
(2.12b) 



On the other hand D A and D$ are the corresponding Green's functions which cannot be separated into two parts by 
removing one tunneling matrix W . Using analytic continuation the current Is can be written as 

I s = -ie|Wf £ J ^ {KH - " K(^) - ^( W )]^( W )} 



dw 



£ / — 2\W\ 2 \D A (u;)\ 2 nt<(u)A s (u;)+lm{t A (cj)}D<(cj) 



(2.13) 



Since 1$ = in steady state, we find the equation relating the lesser part to the retarded part of the Green's function 
or the spectral function As for the dot S. 



Z?<H=7T- 



/eff(w) 



■JmS A {u) 



4s M, 



(2.14a) 
(2.14b) 



That is, the condition that no net current flows out of the dot S determines its nonequilibrium thermal distribution 
function. The current Is can be expressed in another form 



e E / ^ m\ 2 \D A ^)\ 2 7rE<(c)i s H+Im{S^H}^<M 



2tt 



(2.15) 



C. Interacting Case: Ua = and Us 7^ 



When Ua — and Us ^ 0, the active dot A becomes a resonant level and the side-connected dot S is strongly 
correlated. This model Hamiltonian may describe the situation that the dot A lies in the conductance peak region 
and the dot S in the Kondo or Coulomb blockade region. Green's functions of the two dots are given by the Dyson 
equations with the self-energies being given by 



S A (t,t / ) = S c (t, t')+W 2 D s (t,t'), 
X s (t,t') = Xu(t, t') + W 2 D A (t, t'). 



(2.16a) 
(2.16b) 



Here £ c (f, t') — Y, c i,{t,t') + £ c i?(i,i') [equation (2.8a)] is the self-energy coming from tunneling into left and right 
electrodes, Sy(t, t') is the self-energy for a dot S due to the on-site Coulomb interaction, and the auxiliary Green's 
functions are defined by the Dyson equations 

(2.17a) 
(2.17b) 



(2.18a) 
(2.18b) 



D A (t,t') = D A0 (t,t')+ / dti / dh £>A(*,*l)£c(*l,*2)£>A0(t2,t'), 

D s (t,t')=D S o(t,t')+ { dh { dt 2 D s (t,t 1 )Eu{h,t 2 )Dso(t 2 ,t'). 



From the self-energy equations we find 

£<M = 2T L (Lu)f L (tu) + 2T R (Lu)f R (uj) + W 2 D<(lu), 
1]^) - EJ(w) = -2t[T L (uj) + T r (uj)} + \W\ 2 \d s (lu) ~ D a s (uj) 



Inserting all the equations into the expression of the current (2.7), we find 



4 



T p = e E / T {^L{Lo)T R {u)\D^{uj)\ 2 [f p {uj) - f p (co)] + 2\W\ 2 T p (lj)\D a (u>)\ 2 [D<(u>) - 2*f p {u)A S (w)]} (2.19) 
From the condition Is = 0, we have 



D<(W)=7T- 



A s (u) 



2^[r L H/ L H +r fl H/ fl (g;)] + tt\W\ 2 A<(lj) 

r L (Lo)+r R (u)+irW\\2As(u;) 
2Trf eS (u)A s (w), 
r L {u)f L (uj)+T R (Lj)f R (uj) 



x A s (a) 



(2.20a) 



It follows from the relation 1$ = that = 27r/ e fif(w)^4s(o;). The effective Fermi-Dirac function f e g(u>) describes 

t he n onequilibrium thermal distribution at the dot S and also at the active dot A. [It can be shown from the equation 
(2.22).] The current flowing from left to right electrode becomes 



^ I 2tt 

T(w) = AT l (uj)T r {uj)\D a (uj)\ 2 



T{lo) \f R (u - f L (u)], 

tt\W\ 2 A s (lj) 



1 



T l (lu)+T r (lo) 



47tT l (lo)T r (lu) 
T l (lu) + T r {lu) 



A A (to). 



(2.21a) 
(2.21b) 



Note that the current II consists of two contributions, coming from the direct path (L — > A — ¥ R) and the indirect 
path (L^A^S^A^R), which are interfereing each other. 

The current is written down in terms of the transmission coefficient, that is, in Landaucr-Buttiker form. In turn 
the transmission coefficient is proportional to the spectral function of the active dot A. To find the spectral function 
^^(a;) we note that Green's functions of the two dots are related to each other by the equations 



D a (uj) = D\{oj) + \W\ 2 D a (uj)D s (u;)D a (uj) 



(2.22) 



To get the desired spectral function we have to calculate the Green's function D r s at a side-connected dot and the 
auxiliary Green's function D A at the active dot is given by the equation (2.17a). Due to the on-site Coulomb interaction 
at a side-connected dot S, the calculation of D r s {u>) is highly non-trivial and requires a many-body calculation. In 
the next section we describe the non-crossing approximation in order to find the Green's function at a side-connected 
dot S. 



III. NON-CROSSING APPROXIMATION (NCA) 

To study the I-V curves of our model system, we need to calculate the Green's function at a side-connected dot S 
out of equilibrium. For this_pairpose we adopt the K1A which has been a very fruitful tool for the study of Anderson 
model in both equilibriurrJlil'Ej and nonequilibriurrJiSEZI. 

To begin with, we need to identify the effective continuum band for the dot S. In the Anderson model which 
describes the magnetic ions imbedded in normal metals, the self-energy can be written as a sum of two contributions. 

^t. +W (3-D 
J 7r nuj — e + id 

The first term is the one-body contribution from tunneling into the conduction band and the second term E[/ is the 
many-body contribution from the on-site Coulomb interaction. On the other hand the self-energy of the dot S in our 
model is given by 

H s (u;) = \W\ 2 D r A (u) + Tr u (u). (3.2) 



Comparing the two self-energies (equations 3T and ^2) we can identify that the effective Anderson hybridization 
function is given by the equation 
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T eS (e) = -\W\ 2 lmD A (e) 



iWTIm 



= n\W\ 2 A A (e). 



(3.3) 



The effective continuum band for the dot 5* is represented by the auxiliary Green's function D A (t,t'). 

To find the Green's function for the dot S we do a perturbative expansion in W. To begin with we represent the 
Fock space of the dot S by the pseudo particle operators. In the large Us limit we may neglect the double occupation. 



dsa -> b t f a 



(3.4) 



The boson operator b annihilates the empty state and the pseudo fermion operator f a destroys the singly occupied 
state at the dot S. After replacing the electron operator ds a by the pseudo particle operators, the occupation 
constraint has to be enforced: Q = 1. 



(3.5) 



This constraint is taken into account using the Lagrange multiplier method: H — > H + X(Q — 1). The projection to 
the physical Hilbert space Q — 1 leads to the following projection rule for the physical observables: 



0(w) = — lim e f3X O x (uj -> uj + A), 

Zg A— -too 



Z.q = lim e* 3X < Q > 



A — *oo 



A • 



(3.6) 
(3.7) 



Here 0\(lo) is a thermal average of observable operator O for the Hamiltonian H + X(Q — 1). The Feynman diagram 
rules for the Green's functions are very simple. In the 2n-th order there are in + 1 pairs of creation and annihi- 
lation operators. This leads to the factor: [ih]~ 2n [ih] 3n+1 . Canceling one factor ih leads to the following factor: 
( — 1) F ri^]" | W^| 2n where Nf is the number of closed fermion loops. With this rule the NCA self-energies (displayed 
in Fig. 0) are given by the equations 



£/(A;i,t') = ih\W\ 2 G b (X;t,t')D A (t,t'), 
E b (X;t,t') = -ihN s \W\ 2 G f (X;t,t')D A (t',t). 



(3.8a) 
(3.8b) 



Here N s = 2 accounts for the spin degeneracy of singly occupied state for the dot S. Gb and Gj are the ppudo 
particle Green's functions of boson and fermion, respectively. Using analytic continuation to the real-time axisllia we 
find 



X>(X;t,t>) = \W\ 2 G>(X;t,t>)D>(t,t>), 



E<(A;i,i') 

£/(W) 

£>(A;M') 
£<(A;M') 
EE(A;t,f) 



-\W\ 2 G<(X;t,t')D<(t,t'), 
\W\ 2 [Gl(X; t, t')D>{t, t>) - G< (A; *, t')D A (t, t') 
N s \W\ 2 G>(X;t,t')D<(t',t), 
-N s \W\ 2 G<(X;t,t')D>(t',t), 

N S \W\ 2 \G r f (X; t, t')D<(t', t) + G<(A; t, t')D a A (t', t) 



To give a concrete example of analytic continuation we consider the derivation of S^(A; t, t'): 



1 ~ 



-S^(A;i,0 = -ih\W\ 2 -G^(X;t,t') x — D^(t\t). 
i i 1 i 

Removing the common factors on both sides of the equation, we find the desired expression for S^(A; t, t'). 
After Fourier transform, we get the self-energy equations 



(3.9a) 
(3.9b) 

(3.9c) 
(3.9d) 
(3.9e) 
(3.9f) 



(3.10) 



X>(X;lu) = \W\ 2 
E<(A;w) = -\W\ 2 



2n 



G>(A;o;-C)^(C), 



2tt 



G<(A;o;-C)^(C), 



(3.11a) 
(3.11b) 



G 



E>(A;a,)=iV 5 |W| 2 /^G>(A; W + C)fl<(C), 



S<(A;o;) = -iV s |W| 



2vr 



G<(A;c + C)^(C), 



^ b {X-u)=N s \W\ 2 J g [G^(A; W + C)^(C) +G<(A; W + C)^(C)^ • 
The Green's function at the dot S* can be written down in terms of pseudo particle operators. 

ihD s (X;t,t') = <Td Sa (t)4 a (t')> = <Tb\t)f a (t)fl(t')b(t')> . 
In the leading approximation the desired Green's function is 

D S (X; t, t') » ^ G/(A; t, t')G b (X; t', t). 
Analytic continuation to the real-time axisEl leads to 

D>(X;t,t') = -G>(X;t,t')G<(X;t',t), 
D<(X;t,t') = Gf(X;t,t')G>(X;t',t), 

D r s (X; t, t') = ~G r f (X; t, t')G< (A; t' , t) - G< (A; t, t')G a b (X; t' , t). 
After Fourier transform we get 

fdc 



D>(X;lu) 
D<(X;u) 
D r s (X;u) 



2?r 



G>(X;lu + ()G<(X;C), 



^G<(A; W + C)G>(A;C), 



2tt 



G) (A; lu + QG< (A; C) + G< (A; u + C)G£(A; C) 



(3.11c) 
(3. lid) 
(3. lie) 
(3.11f) 

(3.12) 

(3.13) 

(3.14a) 
(3.14b) 
(3.14c) 

(3.15a) 
(3.15b) 
(3.15c) 



From the projection to the physical Hubert space Q = 1 of the Green's function at the dot S, we find the projection 
rules for the pseudo particle Green's functions. 



D^{oj) = — lim e? x D r 8 (\;u) 

ZjQ A— >oo 

= ± l im ^ I ^ 

2tt 



Z s X 
1 



G^(A; w + C + A)G< (A; C + A) + G<(A; to + C + A)G£(A; C + A) 



'GJ( W + 0Gf(0+^(w + 0G?(0 



2vr 

^[iV s G<M-G<H 



(3.16) 
(3.17) 



Since the pseudo particle Green's functions are centered at lu = Es + A (pseudo fermion) and at to — A (slave boson), 
we have to shift the integration variable before taking the limit: A — > oo. For D t s (lu) to be well-defined the limiting 
pseudo particle Green's functions should be defined as 



G r bJ (u)= lim G r bf (X;oj + X), 



G< f (to)= lim e /3A G< / (A;w + A). 



(3.18a) 
(3.18b) 



The second equation means that lim,\— >00 G^(A; lu + A) = 0. Applying the projection procedure to the lesser and 
greater Green's functions of the dot S we find the additional projection rules for the pseudo particle Green's functions: 

Gu Ju) = lim Gg* JX; w + A). (3.19) 



7 



Accordingly the lesser and greater Green's functions of the dot S are given by the equations 

£f H = ^ / § G<(uj + C)G>(C), (3.20a) 

D>(u) = ~j^ G>(co + ()G<((). (3.20b) 

Applying the projection procedure to the physical Green's function D5 at the dot S we obtained the projection rules 
for the pseudo particle Green's functions. 

We are now in a position to find the self-energy equations projected to the physical Hilbert space, Q = 1, for the 
pseudo particle Green's functions. 

Sr^(w) = lim Sr'>(A;o; + A), (3.21a) 
E5 f (w)= lim e /3A E5 f (A;o; + A). (3.21b) 

By comparing the self-energies we can see the identities: Glj(u) = G^j(co). 

Y,<{w) = -\W\ 2 j g G<(cj-()D<(0, (3.22a) 

EJMH^I 2 /^ 0^(0, (3.22b) 

Zf(u) = -N a \W\ 2 J g G<(o; + C)^(C), (3.22c) 

mu;) = N s \W\ 2 J g G5( W + C)^(0- (3.22d) 

These are the NCA self-energy equations out of equilibrium. Now inserting the following equations 

D<(co) = \D r A (w)\ 2 [2T L (u)f L {u) + 2T R (oj)f R (u;)}, (3.23a) 
D>(cu) = \D r A (cu)\ 2 [2T L {u)f L (-u) + 2T R {u)f R (-u)], (3.23b) 

we find the final version of the NCA equations 

S / M = - E / V l^| 2 |^(C)| 2 r p (0 f P (0 Gf{u C), (3.24a) 

S /M = E / - |wTraC)l 2 r P (C) / P (-C) G^(w - 0, (3.24b) 

p=L,R n 

S<H = -iV s ^ |VK| 2 |^(C)|%(0 / P (-C) G<(c + C), (3.24c) 

p=L,R J 77 

SEH = JV a E / v l^l 2 raOI 2 r P (C) / P (C) G r ,(" + 0- (3.24d) 

The effective Anderson hybridization for a side-connected dot S 1 is given by T s(co) = \W\ 2 \D A (aj)\ 2 [T l(uj) + T R (u>)]. 
Although the dot S is not directly coupled to the external electrodes it can experience the bias potential difference 
via the active dot A. 

IV. NUMERICAL RESULTS 

For the numerical calculations, we assume the constant energy-independent tunneling matrix between the active 
quantum dot A and two electrodes. The left and right electrodes are assumed to be described by the same Lorentzian 
density of states (DOS). The self-energies due to the tunneling into the electrodes are found in a closed form in this 
case. The auxiliary Green's function for the active dot A (D' A ) and the effective Anderson hybridization function for 
the side-connected dot S (r e ff(w) = + T r {oj)) are given by the equations 
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D r A (u) = 5—, (4.1) 



hui — Ea — r x 



Tiw + ID 



Here D is the conduction band width of the Lorentzian DOS. The tunneling rate of an active dot A into the electrodes 
r = + r# or the effective band width of Anderson hybridization r e g(w) will be taken as an energy unit in our 
numerical works. For this model the Kondo temperature at zero source-drain bias voltage can be estimated within 
an order of 1 by the equation 

(-fr)- <"> 

Here r e ff = |W| 2 /r is the effective Anderson hybridization at to = and Ea = 0. 

Self- energies: The self- energ y Eg of a side-connected dot S is obtained from the Green's function D$ which is 
calculated from equation (3Tq), using the relation = — [Dg ] _1 . The imaginary part of is displayed in 



Fig ||. Displayed results are corrected to satisfy the causality relation for the Greep^s function at the active dot A. 
This correction procedure is explained in the Appendix C and see also the referencesliS. With decreasing temperature 
a sharp feature develops near u> = 0, which is a sign of Kondo resonance peak. When a finite source-to-drain bias 
voltage is applied the sharp structure near lo = gets smooth and split into two at u> = ±eV/2. 

Spectral functions: The spectral functions are calculated from the imaginary part of Green's functions 

Ai(u) = --ImDJXw), i = A,S, (4.4) 

and displayed in Figs. || and ||. The spectral function As(ui) at a side-connected dot S develops the Kondo resonance 
peak near u) = with decre asing temperature below T < Tk- The Green's function at the active dot A is renormalized 



according to the equation ( |2.22| ). As the side-connected dot S develops the Kondo resonance, the spectral function 
Aa{oj) at the active dot A is depleted near ui = 0. The spectral depletion can be understood as the destructive 
interference between the direct path and the indirect path (Fano interference). The spectral depletion leads to the 
s upp ression of the current flowing through the dot A. The width of depleted region is given by the Kondo temperature 



(4.3) of the dot S. When a finite value of source-to-drain bias voltage is applied, the peak in As(w) and the dip in 
Aa(oS) near oj — get splits into two at ui = ±eV/2. The two Kon do ii re sonance peaks in As(lo) and dips in Aa(u) 
become pinned at the Fermi energies of the left and right electrodes.li^O 

Since the NCA underestimates one-body contribution of the self-energy S5 (side-connected dot) from tunneling 
into the effective continuum band, our NCA calculation of the spectral function at the active dot A shows the negative 
spectral weight near the depleted region (violation of the causality relation). This is a purely artifact of the NCA. 
Since the NCA is based on the expansion in powers of the Anderson hybridization (in our case W) and includes only a 
subset of diagrams up to the infinite order in W without vertex corrections, the underestimation is expected. Though 
the NCA underestimates the one-body contribution of the self-energy, its dependence on the energy and temperatuca 
turns out to be correct close to to — and at T > OK. The NCA leads to the pathological results at and near T = OK.Iiil 
From the Fermi-liquid relation it can be shown that the spectral function of the active dot A is zero at lu = at zero 
temperature. According to the Fermi- liquid theory, the magnetic ion's spectral function (As {to) in our case) can be 
related to the phase shift of conduction electrons (electrons-at the active dot A in our case) at the Fermi energy. The 
phase shift is determined by the Langreth-Friedel sum ruleEj, rig = 2-, where the numerical factor 2 accounts for two 
possible spin directions. 

^(0) = ^-sin^ (4.5) 

7T1 off ^ 

where T e ff = ^S- is the effective Anderson hybridization for a side-con necte d dot S at u = 0. With complete screening 
ns is equal to 1 . Inserting the Fermi-liquid relation into the equation ( 2.22] ) it can be shown that the spectral function 
of the active dot A vanishes at u> = 0, or Aa(0) = 0. 

To compensate the NCA's underestimation of the self-energy we corrected the self-energy by adding a term pro- 
portional to one-body self-energy |W| 2 £)^(o;) to it such that the spectral function Aa of the active dot remains 
non-negative definite and Aa(0) — > with the temperature approaching zero. Since the contribution to the self- 
energy coming from the hybridization with the effective continuum band is of the one-body nature and temperature- 
independent, this correction proedure is not a bad idea. 



9 



Line ar response conductance: The linear response conductance o~(T) is calculated by differentiating the equation 
( 2.21a ) with respect to the source-drain bias voltage at V = 0. 



.(T) = ^/deT(e) . (4.6) 

As shown in Fig. c(T) is suppressed with decreasing temperature. The conductance shows a logarithmic dependence 
on temperature near the Kondo temperature. From Fermi-liquid theory it can be expected that the conductance is 
completely suppressed at zero temperature. The electron flow is blocked by the destructive inteference between two 
different paths: the direct path (L — > A — > R) and the indirect one (L^A^S^A^R). 

This suppressed behavior of a(T) at low temperatures is another example of zero-bias anomalies (ZBA) observed in 
several mesoscopic devices. Scattering of conduction electrons off the Kondo impurities is one of scenarios explaining 
ZBA's. Depending on the topology of mesoscopic devices and Kondo impurites, the ZBA can lead to the enhanced 
or suppressed conductance at low temperature. 

As expected from the Kondo effect, the linear response conductance o~(T) shows the scaling behavior over a wide 
range of temperature. Several curves of cr(T) for different sets of model parameters (differing Kondo temperature) 
collapse onto one curve as shown in Fig. || 



Differential conductance: The current flowing from the left to right electrode is calculated using the equation (2.21a). 
Differential conductance G(T, V) = dl/dV is obtained by differentiating the current with respect to a finite source- 
drain bias voltage. Differential conductance is displayed in Fig. ^. With decreasing temperature dl/dV is suppressed 
near zero source-drain bias voltage. The pseudo gap behavior of the spectral function at the active quantum dot A 
explains this I — V characteristics. 



V. CONCLUSION 



In this work we studied the transport properties of parallel double quantum dots when only the active dot is coupled 
to the two external electrodes and the other dot is side- connected to the active dot. When the active dot lies in the 
conductance peak region and the side-connected dot in the Coulomb blockade region, the current flow through the 
active dot is suppressed due to the Kondo scattering off the effective spin in the side-connected dot. Suppressed 
conductance at low temperature and near zero bias voltage can be understood in terms of Fano interference between 
two paths: direct path of left electrode - active dot - right electrode and indirect path of left electrode - active dot 
- side-coneccted dot - active-dot - right electrode. Our study suggests a possibility that the current flowing through 
one quantum dot can be controlled experimentally by an additional side-connected quantum dot. 
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APPENDIX A: NON-INTERACTING CASE: U A = Us = 



When U a — Us — 0, both active and side-connected dots become resonant levels. The current and its noise can be 
obtained analytically in a closed form. The self-energies of the dots A and S are 

E A (t, t') = E c (t, t 1 ) + \W\ 2 D S0 (t, t'), (Al) 
Z s (t,t') = \W\ 2 D A (t,t'). (A2) 

Here the one-body self-energy S c is given by the equation 

S c (t,t / ) = ^^|V f) (fe)| 2 G' cp (fc;f,t , ) ! 

pk 

and the auxiliary Green's functions D A and D$ are given by the equations 
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D A (t,t') = D A0 (t,t')+ / dh / dt 2 Ai(t,ti)Ec(*i,*a)Ato(*a,*')> 
D s (t,t') = D so (t,t'). 

The self-energy of the auxiliary Green's function Z)^ is given by E^(t, t') = £ c (t, t'). 
The retarded Green's functions can be readily obtained 



(A3) 
(A4) 



D r A (u) 
D r s (cu) 



huj — Es 



[hw - E s ] x [huj -E A - I%(w)] - \W\ 2 ' 

hu - E a -YT c {oj) 
[Huj - E s ] x [Two -E A - E r c (uj)} - \W\ 2 ' 



(A5) 
(A6) 



Due to the Fano interference the spectral function of the active dot A vanishes at u> = Es- On the other hand the 
spectral function of the side-connected dot S is peaked at u> = Es- The nonequilibrium distribution function of the 
dot S is also readily calculated 

E^(w) = 2r i (w)/ i (w) + 2T R (u)f R (u), 



-JmS r A (u)=T Il {u)+r R (u), 



/eff(w) 



r L (u)f L (u) + T R {w)f R {oj) 

r L (w) + r fi ( w ) 



-2ImE^(w) 

The full self-energies of the active dot A are 

£<H = 2T L (u)f L (u) + 2r R (u)f R (u) + \W\ 2 D<(u), 
T, a (lj) - EJ(w) = -2i[r L (w) + r fl (w)] + |VK| 2 - 

To determine Dg(u>), we use the condition Is — (the number of electrons in the dot S remains the same). 

D s {uj) = 2ir As(w). 

r L (w) + t r (uj) 



(A7) 



(A8) 



(A9) 



Here A s (u) = 5(%w - E s ). 

Inserting all the equations into the expression of the current (2.7), we find the current flowing from the left electrode 
to the right one. 

lL=eJ2 ! Trj" {4r £ (w)r fl (w)|I^(w)| 2 [/R(w) - /i(w)] 

a. 

+2\W\ 2 T l (u J )\D a (lu)\ 2 [D<(oj) - 2«f L (fij)A 8 (e>)]} 
= {llj deT(e)[f R (e)-f L (e)]. (A10) 

a 

The transmission spectral function T(e) is given by the equation 



T(e) = 4T L (e)T R (e)\D A (e)f 



7r|W^| 2 i s (e) 
T L (e)+T R (e) 



AnT L {e)T R (e) 
r L (e)+T R (e) 



(All) 



The current noise, defined as the current-current correlation function, can be obtained in a closed form in the nonin- 
teracting case. 



S(t,t') = < 5I(t)5I(t') > + < SI(t')5I(t) > 
Here SI = I— < I >. Defining the current Green's function by 

ihG u (t,t') = < T6I(t)SI{t') >, 



(A12) 
(A13) 
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the current-current correlation function can be written as S(t,t') — hGfj(t,t') + ft,Gfj(t',t). Six diagrams ace con- 
tributing to the current noise and the static current noise Sq = S(u> = 0) is given by the well-known formula.cHI 



So = ^- / d, Tu) !/,„)[! - MmR. /,,>(, i[I -/,,.(<)]} 



4e^ 
' h 



deT(e)[l-T(e)] {f L (e)-f R (e)}\ (A14) 



The first line is the thermal Johnson noise and the second line is called the shot noise coming from a finite source-drain 
bias voltage. 

Typical results at T = OK are displayed in Fig. || for the transmission coefficient T(e) and the Fano factor F = 
So/2eI. Due to the Fano interference the transmission coefficient is completely suppressed at u> = E$, the energy 
level of the side-connected dot. This structure leads to a dip in the Fano factor. The Fano factor F takes the following 
values at V = and V — ► oo. 

F(V = 0) = 1 — ^ 



F(V -> oo) = 1 



(r L + r R ) 2 e% + (e a e s - w 2 ) 2 /r 2 ■ 

2T L T R 



(r L + r fl )2- 

Here T = T L The limiting value of F at eV 3> T is solely determined by the tunneling rates, T L and T R . 



APPENDIX B: INTERACTING CASE: Ua + AND Us + 



When both dots A and S are strongly correlated, the model system becomes highly non-trivial due to many-body 
effects. In this section we derive the effective spin exchange interaction when both dots lie in the Kondo limit or in 
the absence of charge fluctuations. Our model space A4 and its orthogonal space A4 are then given by 



M = \Aa > ®\S/3 > = \Aa;S(3>, a,p = 1=4, 

M = \Aa;Sp >, \Ap;Sa>, \Ap;Sp'>, a, (3 = |>4; P, p' = ±- 

Here a and (5 represent the spin up or down state, while p and p' the isospin state: 

|t+> =4r d lf l 0> ' l*-> = |0 >, i = A,S 

where |0 > means the empty state of a quantum dot. 

To find the effective spin-exchange Hamiltonian we need to evaluate the following matrix elements 

< Aa;S(3\H 1 \Ap;Sp> =0, 

< Aa;SP\H 1 \Ap;Sp> = Sp^ < Aa\H x \Ap > 



pk 



< Aa-Sp\Hi\A P iSp' > = (-lr-'^W 8 a+fjfi 5 p+p , 
When the charge fluctuations can be neglected we find the effective spin-exchange Hamiltonian 



^i = EE Jpp'(k,k') c^ a -a a 

pka p'k' j3 

j pp >(k,k') = 2v p (k)v;,(k') 



i 



-E A E A + U A 



P c P 'k>p ■ Sa 
1 



J AS 



S A - Ss - - 



Jr. 



2\W[ 



1 



1 



E A + U A - E s Es + Us- E A 



(Bl) 
(B2) 



(B3) 

(B4) 
(B5) 



Electrons are scattered off an effective spin S A at the active dot and can flow from one electrode to the other. Since 
S A is coupled antiferromagnetically to Ss, the occurence of enhanced Kondo scattering at low temperature depends 
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on the magnitude of antiferromagnetic coupling J as relative to, e.g., the Kondo temperature. When J as 3> Tk, a 
spin-singlet formation between two dots is favored over the many-body Kondo singlet state and the enhanced Kondo 
scattering is not expected. On the other hand when Tk J as, a Kondo singlet state is expected to be favored 
over a spin-singlet state between two dots and the conductance will be enhanced due to the Kondo scattering at low 
temperature. Detailed study of this model in equilibrium is in progress using the Wilson's numerical renormalization 
and will be published elsewhere. 



APPENDIX C: FERMI-LIQUID RELATIONS 



In this section we discuss the correction of patholog ical b ehavior of NCA near ui = 0. The Green's function of the 



dot A is related to that at the dot S by the equation (2.22). The causality is seldomly violated when the number of 



channels is larger than two. On the other hand the causality is violated for the case of one-channel near u — aX 
temperatures below the Kondo temperature Tk- Recently developed conserving i-matrix approximation (CTMA)EH 
may have a chance to overcome this pathology. In this paper we are going to remedy this pathology of the NCA by 
adding one-body self-energy correction to the self-energy Es(ix>) to satisfy the Fermi liquid relation. The self-energy 
for the dot S can be written as a sum of two contributions 

ESH=Ej yb H + E^H, ShybM = \W\ 2 D r (u;). (CI) 

The self-energy Y7 NCA (bj) calculated from the NCA is known to show the correct energy and temperature dependence 
near lu =r-Q and T = 0. (Temperature must be higher than the pathogical temperature estimated in, e.g., the 
referencestil.) We assume that E[/(w) is well-captured in the NCA. Since — ImEj}(w) ocw 2 + [irksT] 2 for u>,T < 
(look at the referencesE3) , we add the correction term to Ejvca in order to make D r A (ui) obey the causality relation.Ea 

ESM « YT NCA {u) + {\W\ 2 - W 2 NCA )D A (u). (C2) 
Here W^ CA is determined numerically from the calculated Yl r NCA (uj) at the lowest possible temperature. 

= -ImEWO). (C3) 
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FIG. 1. Schematic display of parallel double quantum dots. The active quantum dot A is connected to two left and 
right electrodes and the other dot S is side- connected to the dot A. A side-connected dot S provides an additional current path 
or acts as scattering center of electrons passing through the dot A. 




FIG. 2. NCA self-energy diagrams for (a) pseudo fermion and (b) slave boson. The solid line represents the propagator D, 
the Green's function of the active dot A when W = 0. The dashed (wavy) line is the pseudo fermion (boson) propagator. The 
solid circle means the tunneling matrix W. 
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FIG. 3. Temperature and bias voltage dependence of the self-energy of a side-connected dot S. Temperature varia- 
tion of the self-energy is displayed in (a). The panel (b) is the magnified view of the panel (a) near uj = 0. With 
decreasing temperatures the structure in the self-energy becomes sharper near to — 0. Temperature variations are 
T/T K = 31.6, 10,3.16, 1,3.16 xlO^KT 1 , 3.16 xl(T 2 ,l(r 2 , 3.16 xl(T 3 ,l(n 3 . At T = 1CT 3 xT K the variation of the self-energy 
with increasing bias voltage is displayed in (c). The source-drain bias voltages are varied from eV = to eV = 20 x Tk by the 
amount 2 x Tk. 



1G 






FIG. 4. Temperature and bias voltage dependence of the spectral function As{u>) at a side-connected dot S. (a) and (b) 
show the developments of a Kondo resonance peak with decreasing temperature. The panel (b) is a magnified view of the 
panel (a). The bias voltage dependence of the spectral function As(u>) at T = 10 -3 x Tk is displayed in (c). The Kondo 
resonance peak is progressively suppressed with the increasing bias voltage and develops into two-peak structure. Temperature 
and voltage variations are the same as in Fig. 0. 
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0.3 




FIG. 5. Temperature and bias voltage dependence of the spectral function Aa{u) at the active dot A. Due to the Fano 
interference the spectral function Aa(u>) becomes depleted near lo — as the temperature is lowered below Tk, as shown in 
(a) and (b). When a finite bias voltage is applied, the spectral depletion is refilled and the two-peak structure in As(u>) is 
manifested as the double-dip structure in Aa{w). Temperature and voltage variations are the same as in Fig. ^[ 
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FIG. 6. Temperature dependence of linear response conductance <j(T). The conductance at V = is suppressed with 
decreasing temperature as can be expected from the temperature dependence of the spectral function Aa{co). The conductance 
is suppressed logarithmically due to the enhanced Kondo scattering near the Kondo temperature and becomes saturated with 
decreasing temperature. <j(T) shows the scaling behavior over a wide range of temperature. Symbols mean different sets of 
model parameters or differing Kondo temperatures. 



0.8 - 



CM 

0} 0.6 



CD 



0.4 



0.2 



~\ 1 1 1 1 1 1 1 r 




I i i i i i i i i i I 

01 23456789 10 

eV/T K 

FIG. 7. Differential conductance G(T,V) = dl/dV with varying temperature. The differential conductance is pro- 
gressively suppressed near zero bias voltage with lowering temperature. From top, the temperature variations are 
T/Tk = 10, 3.16, 1, 3.16 x 1CT 1 , 1CT 1 , 3.16 x 1CT 2 , 1CT 2 , 1CT 3 . The last three curves cannot be distinguished with the naked 
eye. 
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FIG. 8. Tunneling coefficients and Fano factor at T = OK. Due to the Fano interference the transmission coefficient 
T(u) is completely suppressed when u> = Es, the energy level of the side-connected dot. The Fano factor, F = So/2eI, 
has a dip structure accordingly. The model parameters are Ea = 0, Es/T — 0.2 for displayed curves. In (a) and (b), 
W/Y — 0.3(1.0) for solid(dashed) line, respectively. In panel (c), the values of W/Y are varied from the bottom successively as 
0.1,0.2,0.3,0.4,0.7,1.0,1.5. 
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